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ABSTRACT 

Ammonium  perchlorate  is  widely  used  as  an  oxidiser  in  composite  propellants  for  rocket 
motors.  Its  decomposition  is  of  interest  in  stability,  sensitivity  and  combustion  studies  for 
rocket  motor  safety  and  performance  analysis  and  prediction.  This  document  presents  some 
theoretical  studies  conducted  to  yield  insights  into  the  mechanisms  and  kinetics  of 
ammordum  perchlorate  decomposition. 
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Ab  Initio  Computational  Chemical  Studies  of 
Molecular  Processes  and  Decomposition  in 
Ammonium  Perchlorate 


Executive  Summary 

The  low-temperature  decomposition  of  solid  ammonium  perchlorate  (AP)  is  a  factor  in 
its  sensitivity,  compatibility  and  safety  and  of  propellants  containing  it.  Furthermore, 
the  products  resulting  from  the  decomposition  can  oxidise  the  binder  and  adversely 
affect  mechanical  and  ballistic  properties.  Consequently,  these  issues  are  of  importance 
to  the  Australian  Defence  Force  (ADF)  for  the  safety  and  service-life  of  munitions. 

In  conjunction  with  experimental  studies,  the  modelling  studies  described  in  this 
report  were  conducted  to  better  understand  the  mechanisms  of  decomposition  in  AP. 
They  will  also  serve  as  a  baseline  for  similar  studies  involving  new  materials  such  as 
ammonium  diniti-amide  (NH4N(N02)2,  ADN)  that  are  proposed  as  oxidisers  in 
propellants  designed  for  low  signature  and  low  sensitivity. 

The  proton  transfer  process  was  studied  in  depth  as  diis  is  thought  to  be  the  preferred 
mechanism  for  the  initial  step  in  AP  decomposition.  In  this  process,  a  proton  (H+)  is 
transferred  from  an  ammonium  cation  to  a  perchlorate  anion  to  form  the  neutral 
species  ammonia  and  perchloric  acid.  For  the  proton  transfer  to  occur,  a  transition 
state  must  be  passed  through  which  is  higher  in  energy  than  either  the  reactants  or  the 
products.  The  energy  required  to  reach  this  state  is  known  as  the  activation  energy. 

Although  AP  is  an  ionic  crystalline  solid,  for  computational  tractabUity,  initial  studies 
modelled  the  transfer  of  a  proton  from  a  single  ammonium  cation  to  a  single 
perchlorate  anion  (this  is  often  referred  to  as  a  "gas  phase"  reaction,  because  it  is  only 
in  the  gas  phase  where  it  is  possible  to  ignore  the  interaction  of  other  ions  or 
molecules).  These  calculations  indicated  that  the  activation  energy  required  for  proton 
transfer  is  very  small  and,  as  might  be  expected,  the  presence  of  ions  in  the  gas  phase  is 
not  favoured  and  perchloric  acid  and  ammonia  is  the  preferred  arrangement. 

Later  calculations  were  conducted  using  a  matrix  of  ions  around  the  ions  involved, 
thus  simulating  proton  transfer  in  the  solid  state.  The  presence  of  the  crystalline  matrix 
was  found  to  inhibit  the  proton  transfer  process  as  the  formation  of  neutral  species 
disrupts  the  bonding  in  the  crystal.  However,  crystal  relaxation  and  the  migration  of 
the  neutral  species  after  proton  migration  would  be  expected  to  lower  the  activation 
energy  of  the  decomposition  reaction,  giving  activation  energies  of  the  same  order  as 
those  observed  experimentally.  Unfortunately,  with  the  computer  resources  currently 
available,  it  was  not  possible  to  model  these  effects.  As  computational  power  increases, 
it  is  planned  to  repeat  these  calculations  using  more  advanced  methods. 
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1.  Introduction 


Conventional  solid  composite  propellants  for  rocket  motors  are  composed  of  particles 
of  an  oxidiser  embedded  in  a  matrix  of  polymeric  binder.  Ammonium  perchlorate  (AP) 
is  widely  used  as  the  oxidiser  material  because  it  is  cheap,  contains  a  large  amount  of 
oxygen  and,  on  combustion,  converts  entirely  into  gaseous  reaction  products.  Some 
properties  of  AP  are  given  in  Table  1  [1]. 

Table  1.  Properties  of  ammonium  perchlorate  [1] 


formula: 

NH4CIO4 

molecular  weight: 

117.5 

density: 

1.95  g/cm^ 

energy  of  formation: 

-2414  kj/kg 

enthalpy  of  formation: 

-2518  kJ/kg 

heat  of  explosion  {H2O  liq.): 

2045  kJ/kg 

oxygen  balance: 

+34.04  % 

nitrogen  content: 

11.04  % 

volume  of  detonation  gases: 

803  L/kg 

decomposition  point  (deflagrates): 

350  °C 

AP  is  a  crystalline  material  composed  of  a  network  of  ammonium  cations  (Nkh^)  and 
perchlorate  anions  (CIO4').  It  is  stable  at  room  temperature,  but  decomposes  at 
measurable  rates  above  150°C  [2]k  In  the  1950s,  the  initial  step  in  the  low-temperature 
(<  300°C)  decomposition  of  AP  was  believed  to  be  one  of  the  following  [3, 4]: 

•  proton  transfer  from  the  ammonium  cation  to  the  perchlorate  anion  to  form 
ammonia  (NH3)  and  perchloric  acid  (HCIO4); 

•  electron  transfer  from  the  perchlorate  anion  to  the  ammonium  cation;  or 

•  Cl-O  bond  breakage. 

Weight-loss  studies  found  that  the  activation  energy  was  the  same  for  sublimation  and 
low  and  high  temperature  decomposition,  being  around  125  kj/mol  [5].  This  lead  to 
the  conclusion  that  a  single,  unified,  mechanism  is  responsible  for  all  these  processes, 
believed  to  be  proton  transfer,  as  shown  in  Figure  1  [5].  However,  if  the  rate  of  reaction 
for  the  high  temperature  decomposition  is  followed  by  pressure  measurements,  an 
activation  energy  of  around  160  kj / mol  is  obtained. 

The  explanation  for  this  follows.  At  both  low  and  high  temperatures,  the  rate¬ 
determining  step  is  the  proton  transfer  that  occurs  in  the  crystalline  AP  to  form 
adsorbed  ammonia  and  perchloric  acid.  At  low  temperatures,  two  adsorbed  HCIO4 
molecules  can  then  react  to  form  chlorine  oxides,  which  then  oxidise  NH3  to  form  the 
ultimate  products  [5,  6,  7].  The  rate  of  decomposition  is  higher  than  that  of  sublimation 


^  However,  newer,  more  accurate  methods,  such  as  microcalorimetry,  may  now  enable  this 
decomposition  to  be  measured  at  much  lower  temperatures,  50°C  or  lower. 
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and  any  gas  phase  reactions  occurring  are  not  rate-determining,  consequently  the  rate 
of  reaction  is  determined  by  the  proton  transfer  with  an  activation  energy  of 
~125  kj/mol,  whether  studied  by  weight  loss  or  pressure  increase. 


high  temperature  mechanism 


Figure  1:  Proton  transfer  mechanism  [5] 


At  higher  temperatures,  the  dissociation  products  are  assumed  to  sublime  before  they 
can  decompose  on  the  surface  [5,  7],  As  the  proton  transfer  is  still  the  rate-determining 
step  for  this  process,  the  activation  energy  measured  by  weight  loss  for  sublimation  is 
the  same  as  that  for  the  low-temperature  decomposition  (~125  kJ/mol).  However,  the 
activation  energy  measured  by  pressure  measurements  gives  the  overall  activation 
energy,  determined  by  the  rate-determining  gas  phase  reactions  (~160  kJ/mol).  In  both 
low-  and  high-temperature  decomposition,  the  ultimate  products  include  O2,  N2,  N2O, 
NO,  CI2,  CIO2,  HCl  and  HNO3  [2,  3,  7] . 

The  low-temperature  decomposition  of  AP  is  a  factor  in  the  sensitivity,  compatibility 
and  safety  of  this  material  and  propellants  containing  it.  Furthermore,  the  products 
resulting  from  the  decomposition  can  oxidise  the  binder  and  adversely  affect 
mechanical  and  ballistic  properties.  Consequently,  these  issues  are  of  importance  to  the 
ADF  for  the  safety  and  service-life  of  munitions.  These  modelling  studies  are  being 
undertaken  in  conjunction  with  experimental  investigations  to  illuminate  the 
mechanisms  of  AP  decomposition.  They  will  also  serve  as  a  baseline  for  similar  studies 
involving  new  materials  such  as  ammonium  dinitramide  (NH4N(N02)2,  ADN),  which 
are  proposed  as  oxidisers  in  propellants  designed  for  low  signature  and  low 
sensitivity. 
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2.  Results  and  Discussion 

2.1  Crystal  structures 

Crystal  structures  were  obtained  from  the  published  crystal  parameters  (described 
below)  using  the  Crystal98  suite  of  programs  [8].  Molecular  visualisations  were 
generated  using  HyperChem®.  The  room-temperature  crystal  structure  of  ammonium 
perchlorate  was  reported  by  Smith  and  Levy  [9]  to  be  Pnma,  a  centrosymmetric  space 
group,  in  accord  with  earlier  studies.  Choi,  et  al.  [10]  also  reported  this  crystal 
structure,  however,  this  was  disputed  by  Peyronel  and  Pignedoli  [11],  who  claimed 
that  the  correct  structure  was  the  non-centrosymmetric  Pnali  space  group.  Choi,  et  al. 
[12]  subsequently  refined  the  structure  of  AP  with  the  X-ray  diffraction  data  of 
Peyronel  and  Pignedoli  but  found  that  the  structure  could  still  be  refined  in  the  space 
group  Pnma.  Lundgen  and  Liminga  [13]  and  Lundgren  [14]  also  found  that  the 
structure  could  be  refined  in  the  Pnma  space  group.  In  any  case,  the  differences 
between  the  two  structures  are  very  small  as  shown  in  Figure  2,  which  shows  the 
crystal  structure  of  ammonium  perchlorate  refined  in  the  PnaZi  space  group  [11]  and  in 
the  Pnma  space  group  at  10  K  [10, 12, 13, 14]. 
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Figure  2:  Crystal  structures  of  ammonium  perchlorate  in  (a)  the  Pnali  space  group  [11]  and  (b) 
the  Pnma  space  group  at  10  K  [10, 12,  13,  14]  (Chlorine  atoms  are  shown  in  yellow, 
nitrogen  -  blue,  oxygen  -  red  and  hydrogen  -  white) 
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The  main  distinction  between  the  Pnali  and  Pnma  space  groups  appears  to  be  slight 
differences  in  the  mean  positions  of  the  ammonium  groups,  which  nevertheless 
undergo  complex  rotational  motions,  even  at  10  K  [9,  10].  As  the  evidence  appears  to 
favour  the  Pnma  structure,  this  will  be  assumed  for  the  remainder  of  this  paper.  In  this 
structure,  each  ammonium  ion  is  surrounded  by  10  oxygen  atoms  with  N-O  distances 
ranging  from  2.9  to  3.5  A  [10],  as  shown  in  Figure  3. 


Figure  3:  Ten  perchlorate  anion  oxygens  (shown  in  green)  nearest  to  the  ammonium  cation  in 
the  Pnma  structure  at  10  K  [10] 

2.2  B3LYP/6-31G(d)  molecular  pair  studies 

Although  AP  is  a  solid  crystalline  material,  it  is  convenient  to  conduct  initial 
calculations  using  an  isolated  pair  of  molecules  or  ions.  These  are  sometimes  called 
"gas  phase"  calculations  as  it  is  only  in  the  gas  phase  where  it  would  be  possible  to 
have  an  isolated  pair  of  molecules  or  ions  interacting.  However,  as  this  is  unrealistic, 
the  effects  of  the  crystalline  matrix  will  be  considered  later. 

Calculations  were  performed  using  a  B3LYP  method  with  a  6-31G(d)  basis  set 
[15, 16, 17]  using  the  Gaussian  98  package  [18]  2.  For  energies  of  optimised  structures,  a 
zero-point  energy  correction  was  applied,  with  a  scaling  factor  of  0.9806  from  Scott  and 
Radom  [20].  This  compares  with  0.9804  given  by  Wong  [21]  and  0.98  given  by 
Bauschlicher  and  Partridge  (0.98)  [22]. 


2  A  general  discussion  on  molecular  orbital  methods  and  basis  sets  may  be  found  in  [19]. 
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2.2.1  Proton  transfer 

2.2.1.1  Potential  energy  scan 

In  the  Pnma  structure  shown  in  Figure  3,  one  hydrogen  atom  of  the  ammonium  ion  is 
pointing  directly  at  one  of  the  oxygen  atoms  of  one  of  the  perchlorate  anions,  as  shown 
in  Figure  4.  This  hydrogen  is  atom  number  9  (hereinafter  called  H9)  using  the 
numbering  scheme  also  shown  in  Figure  4.  It  is  connected  to  the  nitrogen  (atom 
number  6  -  hereinafter  called  N6)  of  the  ammonium  cation,  and  is  pointing  at  an 
oxygen  (04)  of  the  perchlorate  anion.  The  distance  between  H9  and  04  is  1.93  A 
(designated  ro4-H9)  and  the  distance  from  N6  to  04  is  2.95  A,  well  within  the  range  of 
hydrogen  bonding. 


Figure  4:  Ammonium  /perchlorate  combination  in  the  Pnma  crystal  structure  [10]  in  which 
one  ammonium  hydrogen  is  pointing  directly  at  one  perchlorate  oxygen.  Also  shown 
is  the  numbering  scheme  used  to  define  the  atoms. 


This  combination  gives  the  closest  approach  of  a  hydrogen  atom  on  the  ammonium 
cation  to  an  oxygen  atom  on  the  perchlorate  anion  and  appeared  to  be  a  good  structure 
to  start  an  examination  of  the  processes  involved  in  the  initial  phase  of  the 
decomposition  of  AP,  in  which  a  hydrogen  atom  is  transferred  from  an  ammonium 
cation  to  the  perchlorate  anion  to  form  ammonia  and  perchloric  acid.  As  discussed 
above,  this  process  is  thought  to  occur  in  the  AP  crystal  matrix,  however  at  this  stage 
these  species  will  be  considered  in  isolation  from  other  molecules  or  ions,  i.e.. 
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corresponding  to  a  gas  phase  reaction.  The  effect  of  the  ammonium  perchlorate 
crystalline  mati'ix  on  the  reaction  will  be  discussed  in  Section  2.3. 

A  potential  energy  scan  was  undertaken  in  which  the  positions  of  the  atoms  were  held 
constant  in  the  positions  shown  in  Figure  4,  except  for  H9,  which  was  moved 
incrementally  further  from  N6  (increasing  rN6-H9)  and  closer  to  the  perchlorate 
(reducing  ro4-H9).  The  results  are  shown  in  Figure  5,  which  shows  the  energy  relative  to 
the  minimum  energy  with  decreasing  ro4-H9 


Decreasing  r(04-H9)  (A) 

Figure  5:  Variation  in  energy  of  tire  ammonium  /  perchlorate  combination  from  Figure  4  on 
moving  H9  aioay  from  the  ammonium  ion  toioards  the  perchlorate  ion 


It  can  be  seen  that  at  the  high  and  low  extremes  of  ro4-H9,  the  energy  rises  rapidly  as  the 
hydrogen  is  pushed  very  close  to  the  nitrogen  atom  and  very  close  to  the  oxygen  atom 
respectively,  where  coulombic  repulsion  between  nuclei  predominates.  There  are  two 
potential  energy  minima,  one  at  an  ro4-H9  distance  of  about  1.86  A  (rN6-H9  =  1-09  A), 
corresponding  to  the  equilibrium  N-H  bond  distance  in  the  NH4^  /  CIOT  system  and 
another  at  ro4-H9  =  1.06  A  (rN6-H9  =  1-89  A),  corresponding  to  the  hydrogen  atom  having 


3  The  reason  for  defining  ro4-H9  and  not  rN6-H9,  although  not  important  in  this  case,  will  become 
apparent  later  when  a  relaxed  scan  is  performed  on  this  system.  In  a  relaxed  scan,  the  geometry 
is  allowed  to  optimise  at  each  step  (apart  from  the  bond  length  being  defined).  If  rN6-H9  is 
specified  and  ro4-H9  is  allowed  to  optimise,  the  calculations  can  move  the  ammonium  ion  aroxmd 
so  that  H9  (and  rN6-H9  being  restricted)  is  not  necessarily  the  hydrogen  closest  to  the  perchlorate 
ion,  giving  aberrant  results. 


6 


DSTO-TR-1552 


been  transferred  to  the  perchlorate  group  to  give  HCIO4  and  NH3.  In  between,  there  is 
apparently  a  transition  state  (potential  energy  maximum)  at  an  ro4-H9  distance  of  about 
1.43  A  (rN6-H9  =  1.52  A),  giving  an  activation  energy  of  about  38  kj/mol.  As  indicated 
above,  this  is  lower  than  that  expected  for  the  proton  transfer  mechanism 
(-125  kJ/mol)  in  the  crystal.  This  is  not  unexpected,  as  the  proton  transfer  would  be 
expected  to  be  more  facile  in  the  gas  phase,  as  other  ionic  coulombic  attractions  and 
hydrogen  bonds  are  not  disrupted. 

The  results  obtained  from  a  relaxed  potential  energy  scan  on  this  system  are  shown  in 
Figure  6.  In  a  relaxed  scan,  the  geometry  is  allowed  to  optimise  at  each  step  (apart  from 
the  bond  length  being  defined).  Note  that  the  scale  on  the  y-axis  is  much  smaller  than 
that  in  Figure  5.  The  curve  is  very  different  to  that  in  Figure  5,  with  the  overall 
minimum  now  at  an  ro4-H9  distance  of  about  1.06  A,  corresponding  to  the  NH3  /  HCIO4 
system.  The  NITi’^  /  CIO4'  system  has  a  shallow  minimum  at  about  ro4-H9  =  1-68  A  and 
the  activation  energy  is  very  small  (-1.6  kJ/mol). 


Figure  6:  Variation  in  energy  of  the  ammonium  /  perchlorate  combination  from  Figure  4  on 
moving  H9  away  from  the  ammonium  ion  towards  the  perchlorate  ion  allowing  the 
geometry  to  optimise  at  each  step  (except  for  roi-ns  which  is  fixed).  As  discussed  in 
the  next  section,  also  shown  (magenta  circles)  are  the  energies  for  unrestricted 
optimisations  of  the  three  structures  (i.e.  allowing  all  geometric  parameters  to 
optimise,  including  ro4-H9  )■ 
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2.2.1.2  Geometry  optimisation  and  transition  state  structure  identification 

Also  shown  in  Figure  6  (as  magenta  circles)  are  the  energies  obtained  from  unrestricted 
optimisations  on  the  structures  for  the  three  systems,  the  NH4^  /  CIO4'  system  (the 
"initial  structure"),  the  NH3  /  HCIO4  system  (the  "final  structure")  and  the  transition 
state  structure.  The  optimised  structures  obtained  for  these  three  systems  are  shown  in 
Figure  7  and  selected  structural  parameters  are  given  in  Tables  2  and  3. 


(a)  (b)  (c) 


Figure  7;  Optimised  structures  and  numbering  schemes  for  (a)  initial,  (b)  transition  and  (c) 
final  states 
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Table  2.  Selected  angles  between  atoms  ("bond  angles")  of  the  structures  shown  in  Figure  7 


Angle 

Angle  (°) 

Initial  Structure  Transition  State 
Structure 

Final  Structure 

H9-04-C11 

111.0 

108.4 

105.2 

H8-02-C11 

111.3 

113.4 

101.6 

N6-H9-04 

147.3 

172.7 

174.9 

N6-H8-02 

146.2 

118.6 

108.0 

H8-N6-04-C11 

-4.5 

-3.2 

-30.2  4 

H9-04-C11-02 

8.9 

4.6 

46.0 

Table  3.  Selected  interatomic  distances  ("bond  lengths")  of  the  structures  shown  in  Figure  7 


Distance 

— . . . . . B - 

Interatomic  Distance  (A) 

Initial  Structure  Transition  State  Final  Structure 
Structure 

H9-04 

1.66 

1.28 

1.06 

H9-N6 

1.07 

1.23 

1.57 

H8-02 

1.67 

2.17 

2.73 

02-Cll 

1.53 

1.50 

1.48 

03-C11 

1.47 

1.47 

1.47 

04-C11 

1.54 

1.59 

1.64 

05-C11 

1.47 

1.47 

1.46 

H7-N6 

1.02 

1.02 

1.02 

H8-N6 

1.06 

1.03 

1.02 

H10-N6 

1.02 

1.02 

1.02 

N6-C11 

3.23 

3.29 

3.40 

N6-04 

2.62 

2.50 

2.63 

It  can  be  seen  that  the  initial  structure  (Fig.  7a)  is  quite  different  from  that  in  the  Pnma 
crystal  structure  (Fig  4).  In  the  crystal,  the  ammonium  cation  can  hydrogen  bond  with 
several  different  perchlorate  anions  simultaneously;  however,  in  the  initial  structure, 
where  a  single  ammonium  cation  and  a  single  perchlorate  anion  are  isolated,  this  is  not 
possible.  Consequently,  the  ammonium  ion  has  rotated  so  that  two  hydrogen  bonds 
can  form  between  the  ions^.  The  transition  state  structure  is  very  similar  to  the  initial 
state  structure,  accounting  for  their  very  similar  energies  (low  activation  energy).  In  the 
final  structure,  the  NH3  and  HQO4  have  rotated  to  minimise  steric  interaction  and 
maximise  hydrogen  bonding.  This  is  why  this  structure  is  lower  in  energy  than  the 


*  There  is  also  an  alternative  structure  of  the  same  energy  with  dH7-N6-04-cii  =  +30°,  i.e.  where 
atoms  3  and  7  are  involved  in  hydrogen  bonding  instead  of  atoms  2  and  8. 

5  There  is  also  a  stmcture  where  three  hydrogens  of  the  ammonium  are  bonded  to  three  oxygen 
atoms  of  the  perchlorate,  however  this  structure  is  higher  in  energy  than  the  initial  structure 
shown  in  Figure  7a,  probably  because  it  brings  the  nitrogen  and  chlorine  atoms  close  together. 
In  fact,  this  structure  represents  a  transition  state  for  the  ammonium  cation  "tumbling"  aroimd 
the  perchlorate  anion.  This  is  discussed  in  detail  in  Section  2.2.2. 
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initial  structure,  unlike  the  structure  obtained  by  simply  moving  a  hydrogen  from  the 
ammonium  to  the  perchlorate  while  maintaining  the  Pnnin  configuration  as  shown  in 
Figure  5. 

The  initial  structure  has  a  high  degree  of  symmetry,  with  the  angles  aH9-04-cii  and 
aH8-o2-ai  essentially  equal  (Table  2).  In  addition,  the  following  pairs  of  bonds  are 
practically  equal  (Table  3): 

•  the  bond  lengths  between  the  ammonium  hydrogens  and  the  perchlorate 
oxygens  in  the  hydrogen  bonds,  rH9  cm  =  1.66  A  and  rns  02  =  1-67  A; 

•  the  bond  lengths  for  the  N-H  part  of  these  hydrogen  bonds  (i.e.,  the  hydrogens 
pointing  toward  the  perchlorate  anion),  rH8-N6  =  1-06  A  and  rH9-N6  =  1.07  A; 

•  the  N-H  bond  lengths  for  the  H  atoms  pointing  away  from  the  perchlorate 
anion,  rH7-N6  =  rmo-Ne  =  1.02  A.  These  bonds  are,  as  might  be  expected,  slightly 
shorter  than  those  pointing  towards  the  anion; 

•  the  perchlorate  Cl-O  bond  lengths  for  the  oxygen  atoms  facing  the  ammonium 
cation,  ro2-ai  =  1.53  A  and  rcMcii  =  1.54  A;  and 

•  the  Cl-O  bonds  pointing  away  from  the  cation,  ro3-cii  =  roscii  =  1.47  A.  Similarly, 
slightly  shorter  than  those  for  the  oxygens  pointing  towards  the  cation. 

In  the  transition  state  structure,  the  ammonium  group  has  rotated  clockwise  about  an 
axis  perpendicular  to  the  page,  increasing  the  angle  aN6-H904  to  172.7°.  The  dihedral 
angles  dns-NfrCMCii  and  dH9.04.ciio2  are  practically  unchanged,  indicating  that  the 
ammonium  ion  has  not  rotated  relative  to  the  04-C11  bond  and  that  the  H9-04  bond 
has  not  rotated  relative  to  the  rest  of  the  perchlorate  ion.  The  hydrogen  bond 
N6-H8  •  02  has  weakened,  increasing  the  distance  between  H8  and  02  (rH802  =  2.17A) 
and  decreasing  the  distance  between  H8  and  N6  (ms-Ne  =  1.03 A).  On  the  other  hand,  the 
N6-H9  -04  hydrogen  bond  has  strengthened,  decreasing  rH904  to  1.28  A  and 
increasing  rH9-N6  to  1.23 A.  The  bond  lengths  for  the  ammonium  hydrogen  atoms 
directed  away  from  the  perchlorate  anion  (rH7-N6  and  rnio-Ne)  and  for  the  perchlorate 
oxygen  atoms  directed  away  from  the  ammonium  cation  (ros-cii  and  ros-ai)  are 
unchanged. 

In  the  final  structure,  the  hydrogen  has  completely  transferred  to  the  perchlorate  to 
form  perchloric  acid  (rH904  =  1.06 A)  and  ammonia.  The  angle  aN6-H9-04  has  further 
increased  to  174.2°  and  the  dihedral  angle  dH9<M-cii02  has  increased  to  46°  so  that  the 
H9-C)4  bond  is  now  staggered  in  relation  to  the  02  and  03  oxygen  atoms.  The  dihedral 
angle  dns-Ne-CM-cii  has  also  changed  as  H8  has  followed  02  in  its  rotation.  The  distance 
rH8-02  is  too  large  for  significant  bonding  (confirmed  by  the  fact  that  rH7-N6,  rH8-N6  and 
rHio-N6  are  aU  identical),  so  the  attraction  between  H8  and  02  is  probably  almost  purely 
electrostatic  (H8  has  a  net  positive  atomic  charge  of  about  0.36e  while  02  has  a  charge 
of  about  -O.Sle,  see  the  discussion  of  atomic  charges  in  the  following  paragraphs). 
There  is,  however,  stiU  a  hydrogen  bond  N6---H9-04  with  rH9-N6  =  1.57  A,  keeping  the 
pair  of  neutral  molecules  together. 
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Table  4  shows  the  atomic  charges  for  the  atoms  in  the  structures  shown  in  Figure  1,  in 
conjunction  with  those  for  the  atoms  in  the  isolated  species^.  In  the  isolated  ions,  the 
ions  are  totally  independent  with  classical  integer  charge  values  for  the  ions,  -1  for 
CIO4'  and  +1  for  NH4^.  On  forming  the  initial  structure  (the  ionic  complex  NH4^C104’ 
with  two  hydrogen  bonds  N6-H8  ‘02  and  N6-H9-  04),  the  negative  charge  on  the 
perchlorate  anion  decreases  and  the  positive  charge  on  the  ammonium  cation 
decreases  correspondingly  as  electron  density  is  transferred  from  the  perchlorate  to  the 
ammonium.  In  fact  this  is  a  monotonic  trend  from  the  isolated  ions  through  the  initial, 
transition  state  and  final  structures  to  the  isolated  molecules,  when  finally  the  charges 
on  the  two  structures  again  have  classical  integer  values  of  zero. 

In  the  isolated  ions,  the  charge  distribution  is  symmetrically  distributed  among  the 
substituent  atoms,  as  expected.  On  going  to  the  initial  structure,  the  positive  charges  on 
H8  and  H9  have  increased  with  the  formation  of  the  hydrogen  bonds,  at  the  expense  of 
H7  and  HIO.  Similarly,  the  negative  charges  on  02  and  04  have  increased  at  the 
expense  of  03  and  05.  In  the  transition  state  structure,  the  charge  distribution  is 
asymmetric  as  04  and  H9  are  closer  together  than  02  and  H8.  In  the  final  structure, 
after  the  proton  has  been  transferred  and  the  species  rotate  to  a  lower  energy 
configuration,  the  charges  even  out  somewhat  until,  in  the  isolated  molecules,  the 
negative  charge  on  the  perchloric  oxygens  and  the  positive  charge  on  the  ammonia 
hydrogens  are  distributed  fairly  evenly  (apart  from  04  and  H9  which  are  joined,  of 
course). 


Tcible  4.  Atomic  cluirges  for  atoms  in  tire  structures  shown  in  Figure  7  and  for  atoms  in  the 
isolated  species 


Atomic  Charge  (e) 

atom 

Isolated 

Initial 

Transition 

Final 

Isolated 

Ions 

Structure 

State  Structure 

Structure 

Molecules 

totals: 

C104'  -1.0 

C104'  -0.78 

C104'  -0.74 

HCIO4  -0.12 

HCIO4  0.0 

NH4^  1.0 

NH4^  0.78 

NH4^  0.74 

NH3  0.12 

NH3O.O 

Cll 

1.41 

1.47 

1.48 

1.48 

1.49 

02 

-0.60 

-0.63 

-0.59 

-0.51 

-0.45 

03 

-0.60 

-0.50 

-0.48 

-0.48 

-0.45 

05 

-0.60 

-0.49 

-0.48 

-0.46 

-0.43 

04 

-0.60 

-0.63 

-0.67 

-0.66 

-0.60 

H9 

0.46 

0.45 

0.49 

0.50 

0.44 

N6 

-0.84 

-0.90 

-0.92 

-0.95 

-0.89 

H7 

0.46 

0.39 

0.38 

0.35 

0.30 

H8 

0.46 

0.45 

0.41 

0.36 

0.30 

HIO 

0.46 

0.39 

0.38 

0.35 

0.30 

6  Note  the  ordering  of  the  atoms  in  the  table  to  put  04  and  H9  in  the  middle  of  the  table  as  these 
atoms  are  connected  in  the  final  structure. 
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In  summary,  the  proton  transfer  has  resulted  in  the  breaking  of  two  weak  (i.e.,  long) 
hydrogen  bonds  in  the  initial  structure  (H8-02  and  H9-04)  and  the  formation  of  a 
stronger  hydrogen  bond  (H9-N6)  and  a  very  strong  covalent  bond  (H9-04).  Thus,  it  is 
not  surprising  that  the  final  structure  is  lower  in  energy  (albeit  only  slightly)  than  the 
initial  structure.  The  transition  state  represents  an  intermediate  situation  in  which  one 
hydrogen  bond  has  been  broken  (H8-02)  and  the  H9-04  hydrogen  bond  strengthened. 
However,  this  strengthening  of  this  bond  is  not  sufficient  to  entirely  compensate  for  the 
loss  of  the  other,  and  the  transition  state  is  marginally  higher  in  energy  than  the  initial 
state. 

Now  that  we  have  the  three  fuUy  optimised  structures  (i.e.,  stationary  points  on  the 
potential  energy  surface),  we  can  obtain  more  accurate  values  for  their  energies  and 
hence  for  die  activation  energy.  The  process  requires: 

•  running  an  intrinsic  reaction  coordinate  (IRC)  calculation  to  confirm  that  the 
transition  state  is  in  fact  the  correct  transition  state  linking  the  reactants  and  the 
products; 

•  running  frequency  calculations  to  obtain  scaled  zero  point  energies  for  each 
structure;  and 

•  running  high-level  energy  calculations  on  each  of  the  three  structures  to  obtain 
more  accurate  estimates  of  the  energies  involved  including  zero  point  energies. 

2.2.1. 3  Intrinsic  reaction  coordinate  calculation 

An  intrinsic  reaction  coordinate  (IRC)  calculation  was  run  on  the  transition  state  to 
confirm  that  it  connected  the  initial  and  final  structures.  An  IRC  calculation  starts  at  the 
maximum  energy  point  on  the  potential  energy  surface  corresponding  to  the  transition 
state  and  follows  dre  reaction  path  in  both  directions,  optimising  the  geometry  at  each 
point.  In  this  way,  it  indicates  the  reaction  initial  and  final  structures  and  allows 
confirmation  that  the  transition  state  is  the  appropriate  one.  Some  parameters  for  the 
three  structures  discussed  above  are  given  in  Table  6,  along  with  the  final  steps  in  each 
direction  from  the  IRC  calculation. 


Table  5.  Selected  geometric  parameters  for  the  three  structures  shown  in  Figure  7  ivith  those  for 
the  final  steps  in  each  direction  from  tlm  IRC  calcidation  on  the  transition  state 


Parameter 

Initial 

Structure 

Final  step  of 
IRC  towards 
Initial  Structure 

Transition 

State 

Structure 

Final  step  of 
IRC  towards 
Final  Structure 

Final 

Structure 

rH802  (A) 

1.67 

2.03 

2.17 

2.26 

2.73 

rH8-N6  (A) 

1.06 

1.03 

1.03 

1.02 

1.02 

rH9-04  (A) 

1.66 

1.41 

1.28 

1.11 

1.06 

rH9-N6  (A) 

1.07 

1.14 

1.23 

1.43 

1.57 

ro2-cii  (A) 

1.53 

1.51 

1.50 

1.49 

1.48 

ro4-cu  (A) 

1.54 

1.57 

1.59 

1.63 

1.65 

aH9-04-Cll  (°) 

111.0 

109.1 

108.4 

107.1 

105.2 

aN6-H9C)4  (°) 

147.3 

166.8 

172.7 

175.3 

174.9 
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It  can  be  seen  that  in  almost  all  cases,  the  IRC  parameters  are  leading  from  those  of  the 
transition  state  structure  to  the  initial  and  final  structures,  indicating  that  this  transition 
state  is  a  true  transition  state  between  the  other  two  structures.  The  only  exception  is 
the  angle  aN6-H904/  in  which  the  final  step  of  the  IRC  calculation  towards  the  final 
structure  is  not  between  those  of  the  transition  state  structure  and  the  final  structure. 
The  discrepancy  is  minor  however,  and  may  arise  from  inaccuracies  in  the  molecular 
coordinates  from  which  this  parameter  is  derived. 

2.2. 2.4  High  accuracy  energies 

The  energies  of  the  structures  are  given  in  Table  6.  The  total  energy  at  0  K  (Eo)  is  the 
sum  of  the  electronic  energy  and  the  zero  point  energy.  The  total  energy  at  298.15  K  (E) 
is  the  sum  of  the  electronic  and  thermal  energies  (i.e.,  E  =  Eo  +  Evibrationai  +  Erotaaonai  + 
Etransiationai).  The  enthalpy  (H)  and  free  energy  (G)  have  their  standard  definitions,  i.e., 
H  =  E  +  RT  and  G  =  H  -  TS  where  R  is  the  molar  gas  constant  (8.314  kj  mol'^),  T  is  the 
temperature  (298.15  K  in  this  case)  and  S  is  the  entropy.  As  RT  is  the  same  for  aU  the 
structures,  the  relative  enthalpies  of  the  transition  state  and  the  final  structures  are  the 
same  as  their  respective  relative  total  energies  at  298.15  K. 


Table  6.  B3LYP/6-31G(d)  energies  oftlw  three  structures,  "initial  structure",  "transition  state 
structure"  and  "final  structure" 


Energy 

Initial 

Structure 

Energy  (au) 
Transition  State 
Structure 

Final  Structure 

electronic  energy 

-817.876464 

-817.875820 

-817.877704 

zero  point  energy 

0.063665 

0.060211 

0.062361 

total  energy  at  OK 

-817.812799 

-817.815609 

-817.815343 

total  energy  at  298.15K 

-817.805218 

-817.808426 

-817.807428 

enthalpy 

-817.804274 

-817.807482 

-817.806484 

free  energy 

-817.845957 

-817.848545 

-817.849089 

Energy 

Energy  relative  to  Initial  Structure  (kJ/ mol) 

electronic  energy 

0.00 

1.69 

-3.25 

zero  point  energy 

0.00 

-9.07 

-3.42 

total  energy  at  OK 

0.00 

-7.38 

-6.68 

total  energy  at  298.15K 

0.00 

-8.42 

-5.80 

enthalpy 

0.00 

-8.42 

-5.80 

free  energy 

0.00 

-6.79 

-8.22 

As  discussed  above,  in  the  uncorrected  electronic  energies  the  energy  of  the  transition 
state  is  slightly  higher  than  that  of  the  initial  state.  However  once  the  zero  point  energy 
is  added,  the  total  energy  at  0  K  for  the  transition  state  structure  is  apparently  lower 
than  either  the  initial  state  or  the  final  state.  The  author  suspects  that  the  zero  point 
energy  for  the  Transition  State  Structure  may  not  be  accurate,  as  it  is  quite  different  to 
both  that  of  the  Initial  Structure  and  the  Final  Structure,  but,  in  any  case,  the  energy 
differences  between  the  three  structures  are  very  small. 
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22.1.5  Basis  set  superposition  error  (BSSE) 

There  is  one  more  factor  that  still  has  to  be  considered.  The  above  results  are  not 
corrected  for  basis  set  superposition  error  (BSSE).  BSSE  arises  from  the  fact  that  as  the 
two  species  come  closer  together  and  interact  more,  the  electrons  of  each  species  have 
more  of  the  other  species  orbitals  available,  i.e.,  the  basis  set  for  the  transition  state  is 
effectively  larger  than  that  for  the  initial  and  final  states,  leading  to  an  artificial 
lowering  of  the  energy  (see  [19]  for  a  general  discussion  of  basis  sets). 

The  commonly  used  procedure  for  correcting  for  BSSE  is  called  counterpoise 
correction.  This  method  involves  calculating  the  energy  of  each  individual  species  in  a 
complex  calculated  using  all  basis  functions  available  in  the  complex  by  "ghosting"  the 
atoms  of  the  other  species  in  the  complex,  i.e.,  retaining  the  basis  functions  but  setting 
the  atomic  number  to  zero  [18].  The  energy  of  the  complex  is  then  corrected  by 
subtracting  the  difference  between  the  energy  of  the  individual  species  calculated 
using  the  basis  set  of  the  whole  complex  and  that  calculated  using  only  the  basis  set  of 
the  individual  species  [24].  If  we  apply  the  counterpoise  correction,  the  electronic 
energies  relative  to  the  Initial  state  become  2.25  and  -5.69  kj/mol  for  the  Transition 
State  structure  and  Final  Structure  respectively  (compared  to  1.69  and  -3.25  kJ/mol 
from  Table  6).  Consequently,  the  incorporation  of  BSSE  does  not  significantly  affect  the 
above  qualitative  discussion  of  the  results. 

2.2.2  Ammonium  tumbling 

As  mentioned  above,  there  is  a  transition  state  structure  with  three  hydrogen  bonds 
involving  the  hydrogens  on  the  ammonium  cation  and  the  oxygens  on  the  perchlorate 
anion.  This  structure  is  shown  in  Figure  8,  between  two  initial  structures  so  that  the 
type  of  motion  occurring  is  more  readily  apparent.  The  initial  structure  shown  in 
Figure  8(a)  is  identical  to  that  discussed  above  (Fig.  7(a)),  but  has  been  shown  from  a 
different  perspective,  so  as  to  better  illustrate  the  tumbling  process. 

In  the  initial  step  of  the  tumbling  process,  the  ammonium  rotates  about  the  perchlorate 
anion,  leaving  ^e  H9-04  and  H8-02  bonds  intact  and  forming  a  new  bond  between 
H7  and  03  in  the  transition  state  (Fig.  8(b)).  The  ammonium  can  then  either  rotate  back 
to  form  the  original  structure,  or  rotate  while  keeping  the  H9-04  and  H7-03  bonds 
intact  to  form  an  alternative  initial  structure  (Initial  Structure',  as  shown  in  Figure  8(c)). 
Alternatively,  the  ammonium  in  the  transition  state  could  rotate  about  the  perchlorate 
keeping  the  H8-02  and  H7-03  bonds  intact  to  form  another  alternative  initial  structure 
(Initial  Structure",  not  shown).  All  three  initial  structures  thus  formed  would  be 
entirely  indistinguishable,  unless  some  form  of  chemical  labelling  were  used,  such  as 
substituting  one  or  more  of  the  hydrogen  atoms  for  deuterium  atoms  or  isotopically 
labelling  the  oxygens  of  the  perchlorate.  Other  identical  initial  structures  can  be 
imagined  as  the  ammonium  tumbles  about  the  perchlorate.  Selected  geometric 
parameters  for  the  structures  shown  in  Figure  8  are  given  in  Tables  7  and  8. 
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(a)  (b)  (C) 

Figure  8:  Transition  state  structure  (b)  for  ammonium  "tumbling"  between  two  equivalent 
initial  structures  (a  and  c)  and  the  numbering  scheme  used  to  indicate  the 
movement  of  atoms 

The  hydrogen  bond  lengths  rH9-o4  and  rH8-02  increase  on  going  from  the  initial  structure 
to  the  transition  state  structure  as  the  number  of  hydrogen  bonds  increases  from  two  to 
three,  reducing  the  electron  density  available  for  bonding  and  increasing  the  bond 
lengths.  The  bond  length  rH7-o3  decreases,  and  the  bond  lengths  rH7-N6  and  ro3-cii 
increase  as  the  amount  of  hydrogen  bonding  increases  on  going  from  Initial  Structure 
to  Transition  State  Structure  to  Initial  Structure'.  The  other  bond  lengths  do  not  vary  a 
great  deal,  except  that  the  nitrogen  atom  N6  is  closer  to  the  chlorine  atom  Cll  in  the 
transition  state  structure,  as  might  be  expected. 
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Table  7.  Selected  interatomic  distances  (“bond  lengths")  of  the  structures  shown  in  Figure  8 


distance 

- B - 

Interatomic  Distance  (A) 

Initial  Structure  Transition  State  Initial  Structure' 

Structure  for  (estimated  from 

ammonium  Initial  Structure) 

tumbling 

H8-02 

1.67 

1.94 

4.77 

H9-04 

1.66 

2.02 

1.66 

H7-03 

4.77 

1.98 

1.67 

N6-C11 

3.23 

2.92 

3.23 

02-C11 

1.53 

1.52 

,  1.47 

03-C11 

1.47 

1.51 

1.53 

04-C11 

1.54 

1.51 

1.54 

05-C11 

1.47 

1.47 

1.47 

H7-N6 

1.02 

1.04 

1.06 

H8-N6 

1.06 

1.04 

1.02 

H9-N6 

1.07 

1.04 

1.07 

H10-N6 

1.02 

1.02 

1.02 

Table  8.  Selected  angles  between  three  atoms  ("bond  angles")  of  the  structures  shoion  in 
Figure  8 


Angle 

Angle  (°) 

Initial  Structure  Transition  State 
Structure  for 
ammonium 
tumbling 

Initial  Structure' 
(estimated  from 
Initial  Structure 

H9-04-C11 

111.0 

98.3 

111.0 

H8-02-C11 

111.3 

97.8 

46.9 

H7-03-C11 

56.7 

98.0 

111.3 

N6-H9-04 

147.3 

125.5 

147.3 

N6-H7-03 

62.7 

128.2 

146.2 

H8-N6-04-C11 

-4.5 

-51.8 

-112.5 

H9-04-C11-02 

8.9 

57.3 

126.9 

The  charges  on  the  atoms  are  given  in  Table  9.  It  can  be  seen  that  the  major  difference 
in  charge  distribution  on  going  from  Initial  Structure  to  Transition  State  Structure  is  an 
increase  in  both  tibie  negative  charge  on  03  and  the  positive  charge  on  N7  as  they  are 
brought  closer  together.  Similarly,  on  going  from  Transition  State  Structure  to  Initial 
Structure',  both  the  negative  charge  on  02  and  the  positive  charge  on  H8  decrease  as 
the  distance  between  them  increases. 

The  transition  state  for  ammonium  tumbling  has  two  imaginary  frequencies,  indicating 
tliat  it  is  a  second-order  saddle  point  on  the  potential  energy  surface  (i.e.,  a  point  of 
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maximum  energy  for  two  reaction  paths).  The  reaction  paths  involved  can  be  indicated 
qualitatively  by  examination  of  the  normal  modes  of  vibration  for  these  imaginary 
frequencies,  (ire  reaction  path  gives  the  transition  state  structure  for  the  initial 
structure  going  to  the  transition  state  structure  and  back  to  an  initial  state  (ammonium 
tumbling,  as  shown  in  Figure  8).  The  other  appears  to  be  a  continuation  of  the  rotation 
from  the  initial  structure  to  the  transition  state  structure,  bringing  one  of  the  hydrogens 
closer  to  one  of  the  oxygens,  which  may  be  leading  to  proton  transfer  and  would  give  a 
final  structure  similar  to  that  shown  in  Figure  7(c).  However,  as  the  transition  state 
structure  for  ammonium  tumbling  is  higher  in  energy  than  the  one  shown  in 
Figure  7(b),  the  mechanism  discussed  there  for  proton  transfer  would  still  be  preferred. 

Table  9.  Atomic  charges  for  atoms  shown  in  Figure  8 


atom  /  group 

Atomic  Charge  (e) 

Initial  Structure  Transition  State  Initial  Structure' 

Structure  for  (estimated  from 

ammonium  Initial  Structure) 

tumbling 

Cll 

1.47 

1.53 

1.47 

02 

-0.63 

-0.62 

-0.50 

03 

-0.50 

-0.62 

-0.63 

04 

-0.63 

-0.61 

-0.63 

05 

-0.49 

-0.48 

-0.49 

C104  (total) 

-0.78 

-0.80 

-0.78 

N6 

-0.90 

-0.90 

-0.90 

H7 

0.39 

0.43 

0.45 

H8 

0.45 

0.44 

0.39 

H9 

0.45 

0.43 

0.45 

mo 

0.39 

0.40 

0.39 

NH4  (total) 

0.78 

0.80 

0.78 

Table  10.  B3LYP/6-31G(d)  energies  of  the  transition  state  structure  for  ammonium  tumbling 
compared  to  tlw  initial  structure 


Energy 

Energy  of 
Initial 

Structure  (au) 

Energy  of 
Transition  State 
Structure  (au) 

Energy  of  Transition 
State  Structure  Relative 
to  Initial  Structure 
(kj/mol) 

electronic  energy 

-817.876464 

-817.874431 

5.34 

zero  point  energy 

0.063665 

0.064500 

2.19 

total  energy  at  OK 

-817.812799 

-817.809931 

7.53 

total  energy  at  298.15K 

-817.805218 

-817.803955 

3.32 

enthalpy 

-817.804274 

-817.803011 

3.32 

free  energy 

-817.845957 

-817.839720 

16.4 
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2.3  B3LYP/6-31G(d)  studies  of  proton  transfer  in  the  crystalline  matrix 

The  studies  presented  above  consider  an  ammonium  ion  and  a  perchlorate  ion  in 
isolation,  sometimes  referred  to  as  "gas  phase"  studies.  They  indicated  that  there  is 
little  difference  in  energy  between  all  the  structures  considered  and  none  of  the 
activation  energies  obtained  between  initial  and  transition  states  are  close  to  those 
found  experimentally  for  ammonium  perchlorate.  However,  it  is  believed  that  the 
initial  stage  in  the  low  temperature  decomposition  of  ammonium  perchlorate  (i.e.,  the 
proton  transfer  from  the  ammonium  ion  to  the  perchlorate  ion)  occurs  in  the  crystal. 
Here  the  ions  are  surrounded  by  other  ions  and  there  is  likely  to  be  a  much  greater 
energy  penalty  involved  in  disrupting  the  ionic  crystal  by  effectively  replacing  two 
ions  in  the  structure  with  two  neutral  molecules. 

To  study  the  proton  transfer  in  a  crystal  structure,  an  ensemble  of  ions  was  considered 
as  shown  in  Figure  9.  The  proton  transfer  occurs  between  the  two  ions  in  the  centre  of 
the  ensemble,  shown  in  green.  The  surrounding  ions  simulate  the  crystal  environment 
of  the  two  ions  in  the  centre.  This  ensemble  was  considered  to  be  the  minimum  number 
of  ions  that  could  be  used  to  simulate  the  crystal  environment  to  any  degree  of 
accuracy. 


Figure  9:  Ensemble  of  ammonium  ions  and  perchlorate  anions  used  to  simulate  the  environment 
encountered  in  the  crystal  matrix  by  the  centre  ions  (shown  in  green) 
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Unfortunately,  optimising  the  geometry  of  an  ensemble  such  as  this  is  prohibitively 
expensive  in  computer  time,  so  energy  minima  cannot  be  obtained.  Consequently,  a 
scan  sequence  was  used  (similar  to  that  shown  in  Figure  5),  moving  the  proton  from 
the  central  ammonium  ion  towards  the  perchlorate  anion.  The  potential  energy  curve 
obtained  from  such  a  scan  is  shown  in  Figure  10  (shown  with  decreasing  r(O-H)  for 
compatibility  with  Figures  5  and  6). 


Figure  10:  Variation  in  energy  of  the  ammonium  /  perchlorate  combination  in  tlm  crystal 
ensemble  from  Figure  9  on  jnoving  the  hydrogen  of  the  ammonium  closest  to  the 
perchlorate  anion  aivay  from  the  ammonium  ion  toivards  the  perchlorate  ion 


It  can  be  seen  that  as  the  proton  is  progressively  removed  from  the  ammonium  cation 
towards  the  perchlorate  anion,  the  energy  of  the  ensemble  rises.  This  is  to  be  expected, 
because  as  mentioned  above,  the  process  is  creating  two  neutral  molecules,  disrupting 
the  ionic  matrix.  However,  superimposed  on  this  energy  rise  is  a  curve  similar  to  that 
shown  in  Figure  5,  with  two  local  minima,  one  when  r(O-H)  is  between  1.0  and  1.05 
(which  is  also  the  global  minimum  where  the  ions  are  in  their  equilibrium  positions) 
and  one  when  r(O-H)  is  about  1.11  A,  corresponding  to  the  NH3/HCIO4  complex 
formed  in  the  crystal  matrix.  The  energy  required  to  achieve  this  state  is  approximately 
300  kj/mol;  between  these  two  minima,  is  a  local  maximum  corresponding  to  the 
transition  state.  The  activation  energy  required  to  go  from  the  ground  state  to  the 
transition  state  is  approximately  180  kJ/mol,  which  is  of  the  same  order  as  the 
activation  energy  experimentally  observed  for  the  decomposition  of  AP,  which  as 
discussed  in  the  Introduction  above,  is  about  125  kj /  mol. 
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It  is  not  surprising  that  the  experimental  activation  energy  is  lower  than  that  estimated 
here  as  in  the  real  crystal,  relaxation  in  geometry  and  bonding  can  occur  to  compensate 
for  the  proton  ti'ansfer.  Also,  as  mentioned  in  the  Introduction,  the  products  from  the 
proton  transfer  in  the  crystal  are  absorbed  on  the  surface.  Thus,  it  is  likely  that  a 
synergistic  process  takes  place  whereby,  as  the  proton  transfer  occurs,  the  geometry 
and  bonding  in  the  remainder  of  the  crystal  matrix  relaxes  and  the  neutral  species 
migrate  to  tiie  surface.  This  would  have  the  effect  of  lowering  the  local  minimum 
observed  at  r(O-H)  =  1.11  A  (probably  making  it  a  true  minimum)  and  also  reducing 
the  activation  energy.  Unfortunately,  at  this  stage,  computations  to  confirm  these 
suppositions  are  beyond  the  computer  power  currently  available  to  the  author. 


3.  Conclusions 

Aldiough  it  was  not  possible  to  obtain  energies  for  the  various  structures  involved 
which  gave  values  which  were  consistent  with  those  obtained  experimentally,  the 
theoretical  methods  used  did  yield  some  interesting  insights  into  die  possible 
behaviours  of  ammonium  perchlorate  molecules  and  ions. 

In  isolated  molecules  (gas  phase),  there  are  very  small  energy  differences  separating 
the  various  structures  that  may  be  involved  in  proton  transfer  mechanisms  (Initial 
Structure,  the  Transition  State  Structure  and  the  Final  Structure).  As  might  be  expected, 
the  presence  of  ions  in  the  gas  phase  is  not  favoured  and  the  Final  Structure  (i.e., 
perchloric  acid  and  ammonia)  is  the  preferred  arrangement. 

Ammonium  tumbling  around  the  perchlorate  ion  (or  vice  versa)  has  a  very  low 
activation  energy  and  adds  another  avenue  for  atomic  arrangements  that  may  lead  to 
proton  transfer.  Basis  set  superposition  error  (BSSE)  does  not  appear  to  affect  the 
results,  at  least  as  far  as  the  qualitative  discussion  is  concerned. 

The  presence  of  the  crystalline  matrix,  as  would  be  expected,  adds  a  significant  amount 
of  resistance  to  the  proton  transfer  process.  However,  it  was  not  possible  to  model  the 
effects  of  crystal  relaxation  as  neutral  species  migrate  out  of  the  crystal  matrix,  which 
would  be  expected  to  lower  the  activation  energy. 
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Appendix  A:  B3LYP/6-311+G(3df,2p)  Molecular  Pair 

Studies 


Due  to  the  failure  of  the  B3LYP/6-31G(d)  method  to  yield  an  activation  energy  for  the 
transition  state  that  was  consistent  with  experimental  data,  calculations  were 
performed  using  a  larger  basis  set.  Calculations  were  performed  using  a  B3LYP 
method  with  a  6-311+G(3df,2p)  basis  set  [15,  16,  17],  For  energies  of  optimised 
structures,  a  zero-point  energy  correction  was  applied,  with  a  scaling  factor  of  0.989 
[22]. 

A.I.  Rocking  between  two  equivalent  initial  structures 

The  initial  structure  obtained  with  the  B3LYP/ 6-31G(d)  basis  set  (Fig.  7a)  was  used  as 
the  starting  structure  for  a  geometry  optimisation  using  the  larger 
B3LYP/  6-311 +G(3df,2p)  basis  set.  At  this  higher  level,  the  symmetric  structure  shown 
for  the  initial  structure  above  (Fig.  7a)  is  actually  itself  a  transition  state  between  two 
equivalent  structures  in  which  one  N-H  bond  facing  the  perchlorate  anion  is  longer 
than  the  other,  similar  to  the  transition  state  identified  in  Figure  7b,  as  shown  in 
Figure  11. 
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Figure  11:  Transition  state  structure  (centre)  for  "rocking"  between  two  equivalent  initial  state 
structures 


This  may  explain  why  the  zero  point  corrected  energy  for  the  "transition  state"  above 
is  actually  lower  in  energy  than  those  for  the  initial  and  final  states  -  the  transition  state 
is  actually  the  ground  state  for  the  NFl4^  /  CIO4'  system.  Interatomic  angles  and 
distances  are  given  in  Tables  11  and  12  respectively,  using  the  same  atom  numbering  as 
before.  The  energies  for  these  structures  are  given  in  Table  13.  It  can  be  seen  that  the 
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energy  required  to  "rock"  between  the  initial  structures  is  very  small,  around  1.3 
kj/mol  at  298.15  K. 

Table  11.  Selected  angles  between  three  atoms  ("bond  angles")  of  the  structures  shown  in 
Figure  11 


Angle 

Angle  (°) 

Initial  Structure  Transition  State  Initial' 

Structure  Structure 

H9-04-C11 

H8-02-C11 

N6-H9-04 

H8-N6-04-C11 

H9-04-C11-02 

109.6  112.3  114.9 

114.9  113.0  109.6 

168.5  147.6  120.1 

-1.9  -1.9  0.5 

1.5  3.9  -0.7 

Table  12.  Selected  interatomic  distances  ("bond  lengths")  of  the  struchires  shoion  in  Figure  11 


distance 

.  . . . . — n - 

Interatomic  Distance  (A) 

Initial  Structure  Transition  State  Initial'  Structure 
Structure 

H8-02 

2.10 

1.70 

1.41 

H9-04 

1.41 

1.65 

2.10 

N6-C11 

3.26 

3.21 

3.26 

02-C11 

1.46 

1.49 

1.52 

03-C11 

1.43 

1.43 

1.43 

04-C11 

1.52 

1.49 

1.46 

05-C11 

1.43 

1.43 

1.43 

H7-N6 

1.01 

1.02 

1.01 

H8-N6 

1.02 

1.05 

1.13 

H9-N6 

1.13 

1.06 

1.02 

H10-N6 

1.01 

1.02 

1.01 

Table  13.  B3LYP/6-311+G(3df2p)  energies  of  the  new  initial  structure  and  the  structure  for  the 
transition  state  between  tivo  equivalent  initial  structures  as  shown  in  Figure  11 


Energy 

Energy  (au) 

Initial  structure  Structure  for  transition  state 

for  "rocking"  between 
equivalent  initial  states 
(energy  relative  to  Initial 
Structure  (kJ/mol) 

electionic  energy 
zero  point  energy 
total  energy  at  OK 
total  energy  at  298. 15K 
entlralpy 
free  energy 

-818.146785  -818.146416  (0.97) 

0.064164  0.065009  (2.22) 

-818.082621  -818.081407  (3.19) 

-818.075354  -818.074848  (1.33) 

-818.074410  -818.073904  (1.33) 

-818.115277  -818.113189  (5.48) 
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A.2.  Proton  transfer 

A  final  structure  after  proton  transfer  (i.e.,  the  NH3  /  HCIO4  system)  at  the 
B3LYP/ 6-311 +G(3df,2p)  level  can  be  identified  and  appears  similar  to  that  shown  for 
the  lower  level  of  theory  above  (Fig.  7c).  Geometric  parameters  for  this  structure  are 
given  in  Tables  14  and  15.  The  relevant  energies  are  given  in  Table  16. 

The  relevant  geometric  parameters  of  the  NH3/HCIO4  complex  are  compared  with 
those  for  NH3  ^d  HQO4  in  isolation  in  Table  17.  The  parameters  for  HCIO4  are 
essentially  identical  to  those  reported  by  Zhu  and  Lin  [23]  who  studied  the 
unimolecular  decomposition  of  perchloric  acid.  This  is  to  be  expected  as  their 
calculations  were  also  conducted  at  the  B3LYP/  6-311+G(3df,2p)  level. 

Table  14.  Selected  angles  between  three  atoms  ("bond  angles")  for  the  B3LYP/6-311+G(3df  2p) 
final  stnichire  after  proton  transfer  (NH3  /HCl O4) 


Bond  Angle 

Angle  (°) 

H9-04-C11 

107.0 

H8-02-C11 

94.1 

N6-H9-04 

178.7 

H8-N6-04-C11 

-42.3 

Table  15.  Selected  interatomic  distances  and  angles  ("bond  lengths")  for  the  B3LYP/6- 
311+G(3df  2p)  final  structure  after  proton  transfer  (NH3/HCIO4) 


Bond 

Interatomic  distance  (A) 

02-Cll 

1.43 

03-C11 

1.43 

04-C11 

1.60 

05-ai 

1.42 

H7-N6 

1.02 

H8-N6 

1.01 

H9-N6 

1.58 

H10-N6 

1.01 

H8-02 

3.14 

H9-04 

1.05 

N6-ai 

3.44 

The  major  difference  in  the  geometric  parameters  of  the  complex  compared  to  the 
isolated  species  are,  as  might  be  expected,  in  tire  H9-04  and  04-C11  bond  lengths  of 
HCIO4.  In  the  complex,  the  H9-04  bond  is  longer  and  weaker  as  the  hydrogen  is  drawn 
toward  the  NH3  molecule  and,  as  a  consequence,  the  C>4-C11  bond  becomes  stronger 
and  shorter.  Attempts  to  find  a  transition  state  for  proton  transfer  between  the  initial 
structure  for  the  /  CIO4'  system  and  the  final  structure  (the  NH3  /  HQO4  system) 
at  the  B3LYP/6-311+G(3df,2p)  level  failed.  Attempts  using  MP2/6-31G(d)  and 
MP2/6-311G(d,p)  also  faUed. 
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Table  16.  B3LYP/6-311+G(3df,2p)  energies  of  tbe  final  stnichire 


Energy 

Energy  (au)  (energy  relative  to  Initial 
Structure  (kj /  mol) 

electi'onic  energy 

-818.148537  (-4.6) 

zero  point  energy 

0.063726  (-1.1) 

total  energy  at  OK 

-818.084812  (-5.7) 

total  energy  at  298.15K 

-818.077032  (-4.4) 

enthalpy 

-818.076087  (-4.4) 

free  energy 

-818.119062  (-9.9) 

Table  17.  A  comparison  between  geometric  parameters  for  perchloric  acid  in  tlte  NH3/HCIO4 
complex  described  above  with  those  obtained  for  perchloric  acid  in  isolation.  Both 
studies  were  conducted  at  the  B3LYP/6-311+G(3df,  2p)  level. 


Bond  Length  (A) 

NH3/HCIO4 

complex 

NH3  or  HCIO4  in  isolation 

02-Cll 

1.43 

1.42 

03-C11 

1.43 

1.42 

04-C11 

1.60 

1.66 

05-C11 

1.42 

1.41 

H9-04 

1.05 

0.97 

H7-N6 

1.02 

1.01 

H8-N6 

1.01 

1.01 

N10-N6 

1.01 

1.01 

Bond  Angle  (°) 

H9-04-C11 

107.0 

105.9 

04-C11-02 

105.9 

105.1 

05-C11-02 

114.1 

115.1 

H9-04-C11-05 

177.2 

179.7 

H9-04-C11-02 

56.9 

60.4 

H8-N6-H7 

107.7 

107.2 

H7-N6-H10 

108.0 

107.3 

H10-N6-H8 

107.9 

107.3 

H8-N6-H7-H10 

-116.3 

-115.0 
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